ON THE ERROR ESTIMATE FOR CUBATURE ON WIENER SPACE 



THOMAS CASS AND CHRISTIAN LITTERER 

Abstract. It was pointed out in Crisan, Ghazali [2] that the error estimate for the cubature on 
Wiener space algorithm developed in Lyons, Victoir requires an additional assumption on 
the drift. In this note we demonstrate that it is straightforward to adopt the analysis of Kusuoka 
[7] to obtain a general estimate without an additional assumptions on the drift. In the process 
we slightly sharpen the bounds derived in [7|. 



1. Introduction 

In pricing and hedging financial derivatives as well as in assessing the risk inherent in complex 
systems we often have to find approximations to expectations of functionals of solutions to stochastic 
differential equations (SDE). We consider a Stratonovich stochastic differential equation 

d 

(i.i) d&x = Vote,*)* + ° dB i> = x 

i=l 

defined by a family of smooth vector fields Vi and driven by Brownian motion. It is well known 
that computing Pr-tf '■= E(f(^T—t,x)) corresponds to solving a parabolic partial differential equa- 
tion (PDE). The cubature on Wiener space method developed by Lyons, Victoir in IT] , following 
Kusuoka [5] (in the following also referred to as the KLV method) is a high order particle method 
for approximating the weak solution of stochastic differential equations in Stratonovich form. To 
obtain high order error bounds the test functions are assumed to be Lipschitz and the vector fields 
defining the SDE satisfy Kusuoka's UFG condition (see [6]), which is a weaker assumption than the 
usual uniform Hormander condition. 

High order particle methods have since been shown to be highly effective in practice see e.g. 
[15] , [14] and further extensions and applications of cubature on Wiener space have been developed 
by various authors. Applications include the non-linear filtering problem [2], stochastic backward 
differential equations [3], [4], calculating Greeks by cubature methods 17! or extending the KLV 
method by adding recombination in j9] . It was pointed out in Crisan, Ghazali [2] that the analysis of 
the error bounds in Lyons, Victoir [TTJ requires an additional assumption on the drift (see Definition 
I2.5[) and the question was raised if this additional assumption is necessary to derive high order error 
bounds. We first give a brief introduction to cubature on Wiener space and outline how the need 
for an additional assumption on the drift arises in [TTJ . Then, based on Litterer [9] , we demonstrate 
carefully how the analysis in Kusuoka [7] can be adopted to derive similar bounds for cubature on 
Wiener space. We show for the KLV method based on a cubature measure of degree m over a k 
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step partition T> the error E-p can be bounded by 

(2m 
£ ^||/|| V)J . + S (™+1)/2||V/|| 
i=m+l 

for any s € (0,1]. Note that these bounds do not contain any higher order derivatives in the 
direction of the drift Vq and, although our proof contains many elements of the analysis of a version 
of Kusuoka's algorithm carried out in [7], we obtain slightly sharper error bounds in the process 
involving 2m instead of m m+1 derivatives. For suitable families of partitions (first considered in 
Kusuoka [8]) the error bounds immediately lead to convergence of order (m — 1) /2 in the number 
of time steps in the partition. Finally, we clarify the relation of the KLV method to the version of 
Kusuoka's algorithm analysed in [TJ. 

2. CUBATURE MEASURES 

Let C£°(R N ,R N ) denote the smooth bounded R N valued functions whose derivatives of all 
order are bounded. Then V l = (V^,...,Vf) € C^°{R N ,R N ), < i < d may be identified with 
smooth vector fields on R N . Let B — (B}, . . . , Bf) be a Brownian motion and B^(t) = t. Let £t )S , 
t E [0,T], i e R N be a version of the solution of the Stratonovich stochastic differential equation 
(SDE) (11.11) that coincides with the pathwise solution on continuous paths of bounded variation 
(recall that the set of bounded variation paths have zero Wiener measure). We define the Ito 
functional <S> T ,x : Cg([0, T], R d ) -> R N by 

The particular choice for the version of the SDE solution when defining £ t<x implies that the Ito 
functional for a bounded variation path uj coincides with the usual ODE solution of equation (|1.1[) 
along the path u). 

Define the set of all multi-indices A by A = UfcLo{0> ■ • ■ ' d} k an d let a — {ot\, . . . , a>k) € A be 
a multi-index. Furthermore we define a degree on a multi-index a by ||a|| = k + card{j : otj = 0} 
and 

A(j) = {a € A : \\a\\ <j}. 

Let Ai = A\{0, (0)} and A\(j) = {a e A\ : \\a\\ < j}. Following Kusuoka [7J we inductively define 
a family of vector fields indexed by A by taking 

V m =0, Vw = Vi, 0<i<d 

Vf( ai ,...o k ,t)] = [Vj a ],Fi], < i < d,a e A. 
Moreover let V a = V ai ■ ■ ■ V ak where the composition is taken in the sense of differential operators. 
Finally we define a family of semi- norms on the space of functions C£°(R N ) 

k 

ll/lkfc = E £ \\v [ai] ...v [a]] f\u 

j=l ai,...,aj-eAi,||Q!i||H h||«j ||=fc 

It is important to note that these semi- norms contain no derivatives in the direction of Vo- For 
V € C£°(R N ; R N ) we define the flow Exp(tV)(x) to be the solution of the autonomous ODE 

X(t, x) = V(X(t, x)) t > 0, X(0, x) = x e R N . 

A cubature measure on a finite dimensional measure space is a discrete positive measure that 
integrates polynomials up to a certain (finite) degree correctly (i.e. as under Wiener measure). 
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Together with the Taylor approximation for error estimation, cubature is a classical and efficient 
approach to the numerical integration of sufficiently smooth functions. For the Wiener space setting 
[TTj a cubature measure is a discrete measure supported on paths of bounded variation and the role 
of polynomials is taken by the analogous Wiener Junctionals (iterated Stratonovich integrals). 

Definition 2.1. For fixed T > we say that a discrete measure Qt assigning positive weights 
Ai , . . . , A„ to paths 

Wl ,..., Wn eC o V([0,T],i? d ) 
is a cubature measure of degree m, if for all (ij., . . . , ik) G A (m), 

(2.1) e([ odBil---odBil)=J2^ I 

\J 0<ti<—<t k <T J j =1 J0<ti 

where the expectation is taken under Wiener measure. 

By the scaling property of Brownian motion any cubature measure Qt may be obtained from 
Qi by letting uj t -(t) = VTujj (t/T), j = 1, . . . , d and keeping the weights of Q±. 

Taylor expansions play a crucial role in the estimation of the error when we replace the original 
(Wiener) measure by a cubature measure. On Wiener space the bounds for sufficiently smooth 
functions are obtained by considering stochastic Taylor expansion. The following proposition is a 
sharpened version of Proposition 2.1 in 

Lemma 2.2. Let f £ C b °° (R N ) , me N. Then for every t > 

(2.2) /(&,*) = V ai ---V a J{x) ( odB£---odB?* + R m {t,x,f). 

( ai ,..., ak )eA(m) J0< tl <-<t k <t 

And the remainder process R m (t,x, f) satisfies 

rn+2 

sup ^E{R m {t,x,fY) < C V V' 2 sup ||V r ai ...V ai /|| 00) 

x£R N j=m+l (ai,...,c H )eA(j)\A(j-l) 

where C is a constant only depending on d and m. 

Proof. By induction one can prove that the remainder R m (t,x, f) of the Stratonovich stochastic 
Taylor expansion is given by 

R m (t,x,f)= Yl I V ai ...V a J(t tux )odB£...odB%. 

(a, ,„)«(») J0< tl <-<t k <t 

(ui,...,a k )£A(m) 

The proposition follows from an elementary calculation using the Ito formula (see Litterer [TU] for 
details). □ 

The following lemma is the analogous of Proposition 12.21 for the cubature measures Qt and its 
proof may be found in Lyons, Victoir |11) . 

Lemma 2.3. Let R m {T, x, f) be the process defined in 12. 2\) then we have 

m+2 

sup E QT \R m (T,x,f )\ < C(d,m,Qi) V T ] ' 2 sup \\V ai . . . V^/IU, 

x&R N j=m+l (ai,.-.,a z )eA(j)\A(j-l) 

where C is a constant depending only on d, m and the length of the bounded variation paths in the 
support of the cubature measure Q\. 
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The constant in Lemma 12.31 can in fact be made explicit (see Crisan, Ghazali [5] Example 4). 
The expectation of the Taylor approximation /(£t,x) — Rm{s,x,f) defined in (12.21) under Wiener 
and cubature measure coincide by definition of the cubature measure. Hence, one may apply the 
triangle inequality to Lemmas 12.21 and 12.31 and deduce 



h2 

(2.3) sup|£(/(y)-%(/(y)|<c V / 2 sup H^.-.WIIoo. 

xeR N j=m+l (.a u -,on)eA(J)\A(J-i) 

In general, the right hand side of the inequality in (|2.3p is not sufficient to directly obtain a 
good error bound for the approximation of the expectation; in particular, if / is only assumed to 
be Lipschitz the estimate appears useless. Therefore, instead of approximating 

P T f(x) := £(/(&,*)) 

in one step, one considers a partition T> of the interval [0, T] 

t = < h < . . . < t k = T, 

with Sj = tj — tj-i and solves the problem over each of the smaller subintervals by applying the 
cubature method recursively. If r and r' are two path segments we denote their concatenation by 
t ® t' . For the approximation we consider all possible concatenations of cubature paths over the 
subintervals, i.e. all paths of the form uj Si ^ 1 ® . . . ® u Sk ,i k - We define a corresponding probability 
measure v by 

n 

»1 *fc=l 

The iterated cubature method may be interpreted as a Markov operator and, hence the error 
of the approximation of Prf by E v (f(^T,x)) is bounded above by the sum of the errors of the 
approximations over the subintervals. The error over each subintervals can in turn be bounded by 
applying (|2.3[) to Pr-t t f instead of / and exploiting the regularity of Pr-ttf- The following result 
is a corollary to Kusuoka-Stroock [5] and Kusuoka [6 , for a detailed proof see [2] . 

Corollary 2.4. Suppose the family of vector fields Vi, < i < d satisfy the UFG condition. Let 
f € C^(R N ), s € (0, 1] and at, . . . , a, € A\ then 

Cs 1/2 

( 2 - 4 ) II V M ■ • ■ V M P.fU < s( || ai[|+ ;.. +[K[|)/2 HV/||oo 

for all s < 1, where C is a constant independent of s and f . 

As the regularity estimates in the previous corollary do not hold in the Vb direction, but the 
Taylor based estimates used to obtain (|2.3j) require higher derivative in the Vb direction it was 
pointed out in Crisan, Ghazali [5] that the analysis in Lyons, Victoir |IIj requires an additional 
assumption on the drift, we state this as follows. 

Definition 2.5 (V0 condition). A family of vector fields Vi, < i < d satisfies the V0 condition if 

V = u P V lft 

/3GAi(2) 

for some up € C^{R N ). 

The following theorem taken from Lyons, Victoir [TT] is the main error estimate for the iterated 
cubature method. 
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Theorem 2.6. Suppose the vector fields satisfy the UFG and VO conditions then 

i+lfc-i (i+i)/2 



•s 



sup \P T f (x) E v (f(fr, x ))\ <C(T) HV/IU 4 /2 + E E ( T _ , , 

As an immediate Corollary one obtains ([TT] Example 14) high order convergence of the KLV 
method for suitable partitions of [0,T] . 

Corollary 2.7. Consider the family of partitions given by tj = T ( 1 — (l — i) 1 ) and let Vk denote 



3 

the corresponding iterated cubature measures. Suppose the vector fields satisfy the UFG and VO 
conditions then 

sup |Pt/(*)-I^ (/(&,*)) | < Cfc-^-^UV/IU 
where C is a constant independent of k and f. 

In the remainder of the paper we will derive similar bounds for the KLV method that do not 
require the additional VO assumption on the drift. 

3. Algebraic preliminaries - The free Lie algebra and the signature 

In the following we adopt the notation of Lyons, Victoir [TT] . Given a Banach space W we define 
the tensor algebra of non-commutative polynomials over W by 

OO 

T(W) : =0W®'. 

i=0 

Define T^(W) to be the quotient of T(W) by the ideal 0~ J+1 W®\ We identify T^{W) with 
the subspace 

i 

T&(W) = Q)W 9i . 

4=0 

In the following we will not distinguish between the algebras of non-commutative polynomials and 
series as we always work with their truncations. Let eo,...,ed be a fixed orthonormal basis for 
R © R d . Let T(R, R d ) denote the tensor algebra of polynomials over R(B R d endowed with a grad- 
ing that assigns degree two to eo and degree one to the remaining generators (see [TT] for the details 
of the definition). 

Let A e R, a = (ao, Oi, . . .), b — (b , b\, . . .) G T(R, R d ). Define a homogeneous scaling operation 

by 

(A, a) := (a , Aoi, . . . , X l ai, . . .). 

and the exponential and logarithm on T(i?, R d ) using the usual power series. Let ttj denote the 
natural projection of T(R, R d ) onto the subspace T^(R, R d ). 

We define a Lie bracket on T(R, R d ) by [a, b] = a <Eib — b<£) a. Let C denote the free Lie algebra 
generated by R © R d (see Reutenauer [Tj5]). Then C is the space of linear combinations of finite 
sequences of Lie brackets of elements in W — R © R d , i.e. 

W © [W, W) © [W, [W, W}} © • • • . 



6 



THOMAS CASS AND CHRISTIAN LITTERER 



We call an element u of Ttj(C) a Lie polynomial of degree j and an infinite sequence of Lie brackets 
a Lie series. Note that 7Ty(£) C T^(R,R d ). 

Words of the form e a :— e ai <Ei- ■ -<E>e ak , ol G v4u{0} form a basis for T(R, R d ) (note that t% := 1). 
For Wi — X)qgA w ia£a G T(R, R d ), i = 1, 2 we define following Kusuoka [7] an inner product and a 
norm ||-|| 2 on T(R, R d ) by 

(3.1) (wi,w 2 ) = ^ w ^w 2a \\wih = (wi, wi) 1/2 . 

a£iU{0} 

Note that restricted to T^\R,R d ) all norms are equivalent as T^\R,R d ) is finite dimensional 
when regarded as a vector space. 

The map sending to i = 0,...,<i extends to a unique linear map on W and by the 
universality property of the tensor algebra extends to a unique homomorphism F from T(R, R ) 
into the differential operators on The restriction of F to C is a Lie map from C into the smooth 
vector fields on R N . 

Finally we collect a number of simple algebraic facts. 

Lemma 3.1. Let w G C then 

(i) The homogeneous scaling (t,) commutes with exp and log 

(ii) ir m log(ir m w) = ir m log(w) and , n rn exp{'K m w) = -K m exp(w) 

(Hi) r restricted to n m C is a linear map of finite dimensional vector spaces and hence commutes 
with expectations on Tr m £. 

Proof, (i) is obvious from the definition of log and exp as power series, (ii) follows from the fact 
that for a, be T(R, R d ) 7r m (7r m (a)7r m (6)) = 7r m (a&). □ 

For a path <f> G C® bv ({0 7 T], R d ); s,t G [0,T] we define its signature (also known as Chen series) 
S s , t :Cl bv ({0,T},R d j^T(R d )by 

CO „ 

(3.2) S.,tW=S/ #(ti) O • • • ® #(t fc ), 

fc=0^ s< * 1< "' < * fe< * 

where the summation is to be interpreted as a direct sum. Using Stratonovich iterated integrals 
we may define S s j(°B) the random Stratonovich signature of a Brownian motion (under Wiener 
measure) . 

With these definitions in mind we can restate condition (|2.1[) in the definition of a cubature 
measure as 

71 

(3.3) E(ir m (S ,i(°B))) = 2J Aj7r m (S'o,i(w J '))- 

i=i 

Chen's theorem (see e.g. Lyons, Victoir [TT] ) tells us that Li :— ir m (log(So t i(u)i))) is a Lie 
polynomial. The measure Qc = Y^j—y ^j^Lj satisfies 

(3.4) E(Tr m (S ,i(°B))) = E Qc{dL) (TT m exp(L)). 

Conversely for any Lie polynomials Li there exist continuous bounded variation paths with log- 
signature Li. Moreover if Qc satisfies (13.41) Q will satisfy p.3[) , so the identities (|3.3p and (13.41) arc 
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equivalent. The proof of Chen's theorem can be extended to show that log(S s ,t(°B)) is a (random) 
Lie series, see e.g. Lyons [12]. Such arguments can be used to obtain small time asymptotics of the 
solution of Stratonovich SDEs, see e.g. Ben Arous pQ. 

Motivated by this discussion, and following Lyons and Victoir [llj , we make the following equivalent 
definition for a cubature measure on Wiener space. 

Definition 3.2. Let rn 6 N and Qc = Sj=i ^j^L, with A, > and Li € 7r m (£) for i = 1, . . . ,n. 

We say Qc is a cubature measure on Wiener space if and only if 

E(S$(oB)) = E QcidL) w m eML). 

In the following we will sometimes where no confusion arises drop the reference to the integration 
variable L and write Eq c in place of EQ c ( dL y A cubature measure over a general time interval 
[0, T] may be obtained from Q c by homogeneously rescaling the Lie polynomial in its support and 
leaving the weights unchanged. We have 

E(S$(oB)) = E QAdL) Tr m exp((VT,L)). 

4. Error estimate for the cubature approximation 

In this section we derive our main error estimate and demonstrate that Prf can be approximated 
to high order by a cubature measure and the bounds on the error do not involve any derivative in 
the Vo direction (but only its Lie brackets). 

Theorem 4.1. Let P denote the Wiener measure and Q a degree m cubature measure supported 
on paths of bounded variation. Then 

(2m 
E sJ ' /2 ii/ii^-+s (rn+1)/2 iiv/ii 
j=m+l 

for any s £ (0, 1], / G C£°(R N ). The constant C depends on d,m, Q\, 

HVaMOIloo and E P \\ir k (logS(oB))\\l E Ql \\ir k (logS{oB))\\%, 

1 < k < 2m. 

As an immediate consequence we obtain substituting Proposition 3.2 of jTT] by Theorem l4. li the 
following error estimate for the KLV method that preserves its higher order convergence. 



Corollary 4.2. With the notation of Corollary \2. 7| suppose the vector fields satisfy the UFG con- 
dition then 

(k-l I 2m j/2 
4 /2 + E E f '\ (3 -i)/2 +^" +1)/2 

and 

sup \Prf(x)-E Vk {f^ T ,x))\<^~ {m ~ 1)/2 \\^f\\oo 

x£R N 

To prove the theorem we will adopt the analysis of Kusuoka [7j to the cubature on Wiener space 
setting. In the process we sharpen the estimates slightly allowing us to obtain a bound with at most 
2m derivatives instead of m m+1 (compare Kusuoka [7] Lemma 18). Recall that Id is the identity 
function on R N defined by Id(x) — x. 
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£(/(£i,x)) 



Lemma [4 



SDE 



E(Y (TT m S ,i(°B))f(x)) 



T 



Tensor level 



E[ [Y 7r m (exp(-e ) S ,i(°B))] f(Exp(V )(x)) 



Lemma 14.51 



/ (Exp [Y (ir m log(ex P (-E ) S 0tl (oB)))] (Exp(V )(x))) 



Ef(Exp[Tir m logS ,i(°B)](x) 



Flow level 



Before going into technical details we give an interpretation of the ideas developed in Kusuoka [7] , 
summarised in Figure 1. A stochastic Taylor expansion of f(£ x ,s) can be written as T (7r m (S'o,i (°S))) / (x) , 
i.e. the differential operator obtained from the truncated signature under the map Y acting on / 
at x. As the signature is taking values in the tensor algebra we may call the Taylor approximation 
the tensor level (approximation) . It follows immediately from the definition of a degree m cubature 
measure on Wiener space that the expectation of the degree m Taylor approximation under P and 
Q is identical. Although the actual cubature step (exchanging the measures P and Q) has to take 
place at the tensor level we cannot do it directly as the error bounds would involve higher derivatives 
in the direction of Vo, which we have set out to avoid. Instead we follow Kusuoka [7] and observe 
that the signature may be written as the exponential of the log signature and by interchanging exp 
and r we obtain a new approximation at the level of the flow by / (Exp [r (7r m logS'o.i(o_B)))] (x)) . 
Lemma 14.51 formalises this statement and allows us move between the tensor algebra and the flow 
level. Crucially at the level of flows it suffices to approximate £t,x by £t,x m L norm as the bound 
for the approximation of /(£t,x) is only increased by a factor of ||V/||oo- 

A key observation Kusuoka exploits is that if the Lie polynomial defining the flow does not 
involve a eo component the error bound for moving between flow and tensor level does not involve 
higher Vb derivatives. By using a splitting argument at the level of flows (Lemma l4.4[) he can replace 
the log-signature by a Baker-Campbell-Hausdorff style term that does not involve eo. This allows 
him to move to the tensor level and complete the approximation without using higher Vo derivatives. 

To apply Kusuoka's argument to cubature on Wiener space we will go through this approxima- 
tion process (the full lines in Figure 1) for both the Wiener measure and the cubature measure. By 
using the defining cubature identity (|3.3I) we will then be able to see that the approximations at 
the end of each chain agree and obtain the desired bound. 
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The following two lemmas may be found in Kusuoka [7] (Corollary 15 and 17) and we will state 
them without proof. 

Lemma 4.3. Let m > 1, then there exists C > such that 

E P (/(&,«)) - E P [f(Exp [ rvr m (VS, logS 0>1 (oB)) ] (x)) } | < C S ( ro+1 )/ 2 || V/|U 

/or any x G s £ (0, 1] and / G C^{R N ). 

The second lemma is the splitting argument at the level of flows mentioned in the previous 
discussion. 

Lemma 4.4. Let m > 2 and L®, i = 1,2 denote two £( m ) valued random variables with E[\\TTk(L^)\\ 
oo for any k > 1. TTien /or any m > 1 and p G [1, oo), f/iere is C > swc/i £/ia£ 

Bicp(r7r m <>/5,£ (1) )) (W \rn m (^,L^)} (x) 



-Exp 



(w m log(exp(y/s, L {2) ) exp(y/s, £ (1) ))) (x) 



< C- s (™+l)/2 



LP 



for all s G (0, 1] and x G i?^. 
Note that 

£7arp(T7r m (V5,i (1) )) (^arp [^{VJ.LW)] (re)) 
is the composition of Exp(-) functions. 

The following Lemma bounds the difference between flow and tensor approximation. It improves 
on |7J by considering a different truncation of the Taylor approximation. 

Lemma 4.5. Let w = Y^TLi w i € ff m (£), such that each Wi G (tt, — TTj-i) (£), (i.e. each io, is a 
homogeneous Lie polynomial of degree i) and m > 1 tten 

m 

sup |/(Sxp[T(u;)](ar)) - (r [n m exp (w)} f) (x)\ < £ ||r ((7r 2m - 7r„>^') 



// moreover w satisfies (w, eo) = t/ien £/iere is a constant C > smc/i t/ia£ 

2m 

sup |/(£7xp[T(Vi,t«)](a;)-r(7r Tn eip(( > /i,t«))/)(a:)| <C V s J ' /2 ||/|k, 



/or any s G (0, 1], x G i? 



A' 



Proof. We first proceed as in Proposition 9 of [7] by noting that T(w) is a smooth vector field 
defined on all of R N . Thus for any smooth / we have 



-f(Exp(tT(w))(x)) = ((Tw)f)(Exp(tT(w))(x)). 
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Using this identity iteratively to expand f(exp(tT(w))(x)) in a Taylor expansion one sees that 



f(Ex P (tT(w))(x)) - jr t 1 r [ w h ® • • • ® ] ] /(*) 



j'=o 



E 

\-ij<rn, iiH \-ij>m 

Setting t to one we deduce that 



i+-+V=o 

o i! 



[r (w^ ® • • • <g> u^) /] (exp [sT (w)} (x)) ds. 



f(Exp(T(w))(x)) - T 



f E 

\il + --- + ij=0 J ' 



< 



E 



® • • • <g> U>z 



1 Wfy) / 



/(*) 



^2H hij^Tii iiH Mj!>?7t 

m 

<X)||r((7r ain -7r m )^)/L 

3=1 

Noting that n m exp (w) = Ya^+--4 =o jr^i ® ' ' ' ® w b yields the first claim. The second follows by 
considering (y/s, w) in place of w and noting that under the assumption (w, eo) = 0, the vector field 
Vo does not appear on its own in the composition of the differential operators on the right hand 
side of the last inequality. □ 

The last lemma is obtained by combining arguments from Lyons, Victoir [11] and Kusuoka [7]. 

Lemma 4.6. Let t E (0, 1] and Q t be a cubature measure for Wiener space. Then for any x E R N 

E Qt (/(&,x)) ~ E Ql [f(Exp [ Tir m (Vt,logSoA°B))] (x)) } | < Ct( m+1 ^ 2 \\ V/||oo 

for all f E C£°(R N ), where C is a constant independent oft and f. 

Proof. Let g E C£°{R N ). By Lemma HJ we may write 

TO +2 

(4.1) E Qt \g(£ t , x )-r(iT m exp(logS , t (°B)))g(x)\<C V ^'/ 2 sup IKslU 
Letting w = 7r m log Soj(oB) and applying Lemma I4T51 we see that 

2m 

(4.2) £ Qt | 5 (Sxp[L W ](x) - r(7r m e^H) 9 (s)| <CV i J / 2 sup ||V a s||oo 

jtm a€A(j)\A(3-l) 

Combining (|4. 1|) and (|4.2j) with 5 the identity function we see that 



2/1! 



and the lemma follows. 



^ x -Exp[T7T m logS , t (oB)}(x) <CV( )/2 sup ||V a Jd|| 

^ aeA(j)\A(j-l) 



j=m 



□ 
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We are now ready to prove Theorem 14.11 Our proof is modelled along Kusuoka [7] Lemma 18. 
We will go through a sequence of approximations for the expectation of f(£ s ,x) under the Wiener 
measure and the cubature measure. Finally we show that the approximations at each end agree. 

Proof. (Theorem 14. II) 

Let yUi. s = P be the Wiener measure on paths parametrised over [0, s] and /i2,s — Qs- From 
Lemmas 14.31 and 14.61 we see that for i — 1 , 2 



(4.3) sup E» iB {f{t s , x ))-E^Af(Exp[TK m {yrs^gS ,x{oB))]{x))) < C S (" i+1 >/ 2 || V/|| 



Let 

LW = ir m \ogS { m 1 \oB) and = -e„. 

It is well known (see e.g. Lyons [H]) that the log-signature of Brownian motion is a Lie series with 
probability one. Also we have E[\\TTk(logS ^i{oB))\\2] < oo, in fact using the techniques of rough 
paths a similar statement can be obtained at the level of paths. For example Lyons and Sidorova 
compute in |13j the radius of convergence for the log signature. 



Hence, Lemma |4~41 implies that for i = 1,2 



(4.4) 



E^ j(Exp{Tn m {^~s,L^)) (Exp [r7r m (Vi,i (2) >] (*) 



j(Exp Tir m log (exp(y/s,L^) exp{y/s,L^)^ (z) 
x = Exp (T{y/s, -e )) (z) 



< Cs (m+1)/2 ||V/|| 



Writing 

the inequality (|4.4j) becomes 

(4.5) E^j(Exp{Yir m (^- S ,L^)){x) 

- E^ j(Exp [TTT m log (exp{^, exp{^~s, L«))] {Exp (T{^, e » (a))) 

< Cs (m+1)/2 ||V/|| c 

Thus, combining inequalities (|4.3I) and (I4.5[) and using the triangle inequality we see that 

(4.6) sup E^J(C S , X ) 
xeR N 

- E^ j[Exp \rTr m log(exp((^,-e )) (VS, Sb,i(°B)))] (Exp (T{^,e )) (a;))' 

< Cs (m+1)/2 ||V/|| c 

It follows from the Baker-Campbell-Hausdorff formula that 

(4.7) 7T m log(exp((y/s, -e )) (Vs, S ,i(oB)) 
has no eo component, i.e. 



(TT m log(exp((y/s, -e )) (\/5, Sb,i(°-B)}) > e o 



= 0, 
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where we defined the inner product in fl3.1[) . 

Moreover as the log signature of the Brownian motion is a Lie series with probability one, (|4.7p is 
a Lie polynomial. 

Hence, we may apply Lemma 14751 to inequality (|4.6j) and once again using the triangle inequality 
we obtain for i = 1, 2 

(4.8) sup E^ i3 (/(£ s , x )) - -E/w.i ((r7r m exp( v^, 7r m Zo5(exp(-e )5 ,i(oB)) 



<C E ^ /2 ||/lki + s (m+1)/2 ||V/|| 




where r/ = Exp (r(y / s, en)) (z). Note that in the previous step we have used the fact that the scaling 
operation (s, •) commutes with log and exp. We have also used the fact that 

7r m exp ( y/s, w m log(exp(-e ) 5 ,i(oB))) = 7r m (Vs, exp(-ea)So,i(°B)). 

Finally using the cubature relation (|3.3j) 

-Bp (7r m 5*0,1(05)) = E Ql {-K m S ,i(oB))) 

and noting that the multiplication by a deterministic tensor can be taken out of the expectation 
we have 

E P (ir m (exp(-E )S ,i(oB))j = Eq 1 (ir m (exp(—eo)Sa t i(oB))\ . 
Hence, it follows that 

E P [(rir m (^,exp(-e )S 0A (oB))f){y)] = E Ql [(Tn m (yfi, exp{-e )So,i(°B))f)(y)] . 

Using this identity in ()4.8[) a final application of the triangle inequality completes the proof of the 
theorem. □ 

Remark 4.7. The truncated log signatures of the cubature paths of a degree m cubature measure 
satisfy the definition of a m-C-moment similar random variable of Kusuoka 7 with respect to the 
truncated log signature of the Brownian motion. Conversely for any finite such family we can find 
paths that satisfy a degree m cubature formula. The approximation operator for P s f whose error 
bounds are analysed in [7] can be written as 

E Qcf (Exp[T(^/s 7 TT m L)](x)) 

and it is clear from our discussion that the same bounds as in Theorem \4-l\ can be obtained for this 
approximation. 
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